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I describe the computation of energy widths of nuclear states using an integral over the interaction 
region of ab initio variational Monte Carlo wave functions, and I present calculated widths for many 
states. I begin by presenting relations that connect certain short-range integrals to widths. I then 
present predicted widths for 5 < A < 9 nuclei, and I compare them against measured widths. They 
match the data more closely and with less ambiguity than estimates based on spectroscopic factors. 
I consider the consequences of my results for identification of observed states in 8 B, 9 He, and 9 Li. 
I also examine failures of the method and conclude that they generally involve broad states and 
, variational wave functions that are not strongly peaked in the interaction region. After examining 

■ bound-state overlap functions computed from a similar integral relation, I conclude that overlap 

' calculations can diagnose cases in which computed widths should not be trusted. 
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I. INTRODUCTION 



i 1 The last decade and a half have seen enormous progress in the description of light nuclei as collections of interacting 

, nucleons with the same properties as in vacuum [1-4]. After the formulation of potentials that fit nucleon-nucleon 
scattering data with high accuracy, ab initio calculations of nuclear structure demonstrated that the energy spectra 
q \ of nuclei small enough for converged calculations (mass number A < 12) can be understood as arising from the 
vacuum nucleon-nucleon interaction [5j . The addition of a three- nucleon interaction with only two to four parameters 
produces quantitative agreement with experiment for the nuclear binding energies and spin-orbit splittings, including 
the correct ground state (J™ = 3 + instead of 1 + ) for 10 B [6[. 
^-H ' For nuclei with A < 4, methods using correlated hyperspherical harmonic bases or Faddeev and similar formulations 
. have solved bound-state, scattering, and reaction problems quite successfully 0413] ■ The ab initio methods that have 
^sO ' been developed for A > 4 nuclei are suited mainly to treatment of bound states, but there has been significant progress 
on unbound states in recent years. The Green's function Monte Carlo (GFMC) method has been formulated as a 
particle-in-a-box method and used to compute phase shifts in the 5 He system with a full realistic Hamiltonian . 
The no-core shell model (NCSM) has been merged with the resonating group method (RGM) to produce phase shifts 
and reaction cross sections in several systems using effective two-body forces without explicit three-body terms [l~2| — 
[Tf3 |. The coupled-cluster method has been combined with a Gamow shell model basis to compute complex energies 
E = Eft — iT/2 of resonances and thus their excitation energies Er and total widths T. It has been applied to 
resonances of He isotopes and A = 17 nuclei, but it has not been used to produce cross sections [HI, [lTj- All of these 
approaches to scattering are computationally challenging, and significant human and computer effort is required for 
each individual system. 

Even before the development of true scattering/reaction calculations in A > 4 nuclei, there were accurate energy 
calculations for many unbound resonance states [H, H|, produced by approximating them with bound (i.e., square- 
integrable) wave functions and real energies. This approach is successful for states sufficiently narrow that important 
features of their structure can be accommodated within the model space. One way to understand these calculations 
is as approximations to Gamow's decaying complex-energy states [181 ] with energies near the real axis (i.e., small 
widths). 

Quantum Monte Carlo (QMC) calculations of square-integrable "pseudobound" states begin with a variational 
Monte Carlo (VMC) calculation, in which a complicated but closed-form wave function containing large amounts of 
correlation is produced by minimizing the energy expectation value with respect to many variational parameters [l[. 
The computational effort in the VMC method lies in computing energy expectation values by Monte Carlo integration 
over all particle coordinates, and a square-integrable wave function is necessary for both energy and normalization 
integrals to converge. The variational ansatz incorporates square integrability through particle correlations that decay 
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exponentially at large separation, and reasonable energies result because the long-range tails of the true resonant wave 
function have small amplitudes relative to the "interaction region" where all nucleons interact. 

The second step in application of the QMC method to nuclei is a GFMC calculation. This step also requires square- 
intcgrable wave functions for evaluation of matrix elements. GFMC takes the VMC wave function as a starting point 
and evolves the Schrodinger equation for imaginary values of the time t = it through a series of small steps. This 
evolution filters high-energy contamination out of the wave function and for large r leaves behind the lowest-energy 
eigenstate contained in the VMC starting point. For bound or narrow states, this procedure "converges" at large 
t to a unique energy [f| [lj| [2(| . For broad states the GFMC propagation drifts slowly toward the lowest available 
threshold; even in very long GFMC propagations, the curve of energy versus r for such states fails to flatten out at 
a "converged" energy. Presumably this also occurs for the narrow states, but with a decline too slow to have been 
noticed in calculations so far. 

Because the calculated energies of narrow resonances in the QMC methods are believed to be accurate, it is 
natural to ask whether widths could be extracted from these calculations, avoiding the considerable effort of explicit 
scattering/reaction calculations. Attempts to correlate the width with the rate of decline of the energy, dE/dr, in the 
late-time GFMC propagation failed, so the rate may depend more on the starting wave function than on the physical 
width [Hj]. 

Even in explicit scattering calculations, some widths are effectively too narrow to resolve by QMC methods. Direct 



so 



QMC calculations of widths so far consist of energy calculations for a series of specified boundary conditions 11 
the smallest width that can be computed must be larger than the energy resolution of the method. As an extreme 
example, the ground state of 9 B has a total energy of —56.3 MeV and a width of 0.5 keV, roughly 0.01% of the 
total. Since the typical precision claimed for GFMC energy calculations is 1%, this width will not be amenable to 
calculation in this way. 

Here I present a method to extract approximate widths from QMC wave functions. The basis of this approach 
is a relation between the partial width in a specified breakup channel and an integral over the "interaction region" 
where all of the nucleons are close to each other. This relation has been used in the literature, though not (to my 
knowledge) in explicit application to many-body wave functions. It has been used particularly in models of alpha 
[22l I23] ] and proton [24H29 ] decays of heavy nuclei in order to avoid direct integration of the Schrodinger equation 
to impractically large radii. It is also closely related to Green's function formulations of scattering theory. There 
has been considerable interest in applying such formulations to ah initio calculations in the recent literature: they 
have been used to extract asymptotic normalization coefficients (ANCs) from many-body bound states computed in 
various approximations [30Tj33| and are coming into wider use for scattering problems [&, ;34-38|. Preparation for 
future reaction calculations of that kind is a primary motivation of the present work. 

Although the integral relation can be used to compute the partial width in any two-body decay channel if wave 
functions of the parent and daughter nuclei are known, I confine my consideration in the present study to nucleon 
emission. I compute widths for all narrow (T < 1 MeV) one-nucleon decays to bound states in 5 < A < 9 nuclei. I 
also present results for several specific broad states for which better theoretical information concerning widths would 
be useful and for some with unbound final states that are well-represented by pseudobound wave functions. 

The integral relation applied here can be used with either GFMC or VMC wave functions, extending even to use of 
the same computer routines. Application to GFMC requires more computation and additional bookkeeping, so I have 
chosen in this initial study to use only VMC wave functions. This work represents the extension to unbound states of 
the methods presented in Ref. [33| . The integral relation used here may be used with other many-body methods, but 
it is particularly well-suited to QMC wave functions because of its formulation as a short-range integral and because 
of the transparent treatment of fermion antisymmetry in the QMC methods. 

The remainder of this paper is organized as follows: In Sec. [TT] I motivate and define the integral method, and I 
discuss its accuracy and its connection to overlap calculations. In Sec. IIIII I describe the application of the integral 
method to VMC wave functions. In Sec. IIVI I present the results of calculations for specific states, compare computed 
and experimental widths, and consider how accurately widths may be computed from spectroscopic factors. Along 
the way I present new (and likely improved) overlap calculations of VMC wave functions that may prove useful for 
treatment of transfer reactions. Finally in Sec. [V] I summarize my results and briefly mention future directions for 
this work. 

II. AN INTEGRAL RELATION FOR RESONANCE WIDTHS 
A. The connection between widths and asymptotic normalizations 



Consider a many-body wave function "J at energy E above the threshold for breakup into clusters 1 and 2 that 
have wave functions ipi and %p2 respectively and no internal angular momentum. Assume also for the moment that 
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only this breakup channel is open. Since ipi and ip2 are spinless, the orbital angular momentum I of their motion 
is equal to the total angular momentum J of ^. Given appropriate boundary conditions, at large separations 
r\2 of the clusters is a linear combination of an incoming wave proportional to Ii(rj, kr 12) /t 12 and an outgoing wave 
proportional to Oi(j],kri2)/ri2, each normalized to probability flux Hk/p in its appropriate direction at 7*12 — > 00. 
(For all special functions, I follow the conventions of Refs. [39| and [4CJ.) These functions depend on the energy E 
through the wavenumber k = ^/2jlE/h, where p, is the reduced mass of the clusters, and through the Sommerfeld 
parameter rj = ZiZ 2 e 2 p/hk, where Z\ and Z2 are the charge numbers of the clusters. Oi and solve the radial 
Coulomb-Schrodinger equation 

d 2 ui { 1(1 + 1) 2 V \ 

with p = kr\2 and outward and inward flux, respectively, at t\i — > 00. In terms of /; and O;, ^ (assumed to have 
angular-momentum projection to, omitted from subsequent expressions) is, 

#(ri2 00) = C^x^Yimi^) [Ii(v,kr 12 ) - Si(k)Oi(r), kr 12 )] /r 12 . (2) 

where Y/ m (fi 2 ) is a spherical harmonic. In this single-channel case with both incoming and outgoing waves, the 
normalization Ci is arbitrary. If k and 77 are real, conservation of probability guarantees that the effect of scattering 
is to multiply the outgoing wave O; by a complex phase factor relative to the incoming wave Ii, so that the function 
Si(k) may be written in terms of a real phase shift Si(k) as Si — e 2lS ' . The function Si(k) is the single-channel case of 
the S'-matrix, which gives the amplitude and phase of outgoing waves in terms of specified ingoing amplitudes, and it 
determines the cross section uniquely. 

Resonances in the scattering of clusters 1 and 2 occur at real energies near poles of Si(k) and produce peaks in the 
scattering cross section where tan^j ~ T/(E — En) around some resonance with energy Er and width V. S-matrix 
poles in general occur at complex k (and thus complex E). It is evident from Eq. ([2} that at such a pole the wave 
function takes the form 



$(ria -> 00) = C i ViV'2yi( r i2)Oi(TOfcri 2 )/ri 2 , (3) 

with only outgoing flux. This is the same sort of decaying complex-energy state originally formulated in Gamow's 
treatment of alpha decay (l8| . The probability flux out of this state is Hk\C[\ 2 / p,, so the normalization constant is 
no longer arbitrary but carries information about the size of the outgoing flux relative to the total wave function. 
If ^ has been normalized to unit total probability in some finite region, its lifetime is inversely proportional to the 
outward probability flux. By writing ^ in the form 

¥(ria -*■ 00) = CfafoMhz) [[Silk)]- 1 h - O z ] /r 12 , (4) 

and defining to have unit norm in some regularization scheme, it can be shown that the residue of the S'-matrix pol e 
is proportional both to the squared normalization constant \C[\ 2 and to the imaginary part of the pole energy |4dl443| . 
The relations among |C/| 2 , the lifetime, and the pole location imply that the width T = h 2 k\C[\ 2 / p. (Derivations that 
deal rigorously with the normalization of ^ may be found in Refs. [4ll l42j.) The width may thus be computed from 
the normalization constant C[. 

Physically realized systems have real energies, so formulation of the width in terms of complex-energy Gamow states 
is often inconvenient. For the QMC methods, not only is the energy real, but the wave functions are stationary waves, 
sums of incoming and outgoing waves with zero total flux in each channel. Again assuming a single open channel of 
given angular momentum, a standing-wave solution is asymptotically 

*(r X2 -> 00) = a^Ytivu) [Fifakru) + KiityGifakru)} /r 12 , (5) 

for some different Ci than before. Fi(r), kr 12) is the regular Coulomb function that satisfies ^(77,0) = 0, Gi(r],p) 
satisfies the Wronskian relation 

dFi(ri,p) dGi(rj,p) 

— -3 Gi{rj, p) - Fi(r), p) = 1 , 

dp dp 

and the two are related to Oi and by 



Oi(v,p) = G l (r),p)+iF l (r ] ,p) 
Ii(r),p) = Gi(r),p)-iF l (r],p). 



(6) 
(7) 
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The function Ki(k) is the single-channel case of the if- matrix. Eqs. ([2]), ([6]), and (J7J combine to give 

so that Ki(k) = tan Si(k). Eq. ([5]) may then be written as 

*(n 2 oo) = C/Vi^^^ia) [cos^F z (?7, fcri 2 ) + smSiGi(i], kr 12 )} /r 12 (9) 

for some C[ determined by an appropriate normalization of vp. At a resonance, <5; = 7r/2 and dSi/dE > 0. This 
corresponds to a pole of Ki(k) on the positive real axis, so that 

*(r 12 ^oo) = C[M 2 Yi{v 12 )Gi(r,,kri2)/r 12 (10) 
= ^C i ViV'2^(f 12 )[7 / (r ? ,fcr 12 ) + O i (??,A:ri2)]/ri2. (11) 

The standing-wave solution at resonance is thus the sum of equal inward and outward fluxes, each of magnitude 
hk\Cl\ 2 This is the rate at which ^ decays through the outgoing wave and is replenished by the ingoing wave. 
In analogy with the complex- energy case, it can be shown that the residue of the if -matrix at the resonance pole 
is proportional to \C[\ 2 for a wave function normalized to unity in a finite region [44T - |46j . up to corrections arising 
from the choice of normalization volume. These corrections are small as l ong as the wave function has a much smaller 
amplitude in the asymptotic region than in the interaction region [42, 44, 45], as is expected for a long-lived resonance 
state. Comparing Eqs. ([3]) and (|10|). it is unsurprising that the result of a rigorous derivation is T ss h 2 k\C[\ 2 / 'fi. 

I have sketched the considerations that lead to the connection between resonance width and wave function asymp- 
totic normalization for the case of a single open channel and spinless daughter nuclei. These considerations carry over 
directly to ab initio wave functions, but it is necessary to account for complications absent from the toy model. 

Most complications of the ab initio case amount to additional bookkeeping implied by multiple final-state channels, 
antisymmetry of the wave function, and daughters with nonzero angular momenta. Multiple types of decay products 
may be emitted, and their non-zero internal angular momenta may couple in multiple ways. The simple right-hand 
sides of Eqs. @ and <j5j> are replaced by sums over all breakup and I channels of a given total angular momentum J 
and parity tt. The formalism to describe multi-channel wave functions in this case may be found in many treatments 
of reaction theory (e.g. Refs. [13, IH, S3, El] ) • The coefficient Si(k) in Eq. ^ is replaced by a matrix S a b that gives 
the outgoing flux in channel a produced by unit incoming flux in channel 6, and Ki(k) in Eq. ([5]) is similarly replaced 
by a matrix K a b that gives the irregular-function amplitude in channel a produced by unit regular-function amplitude 
in channel b. If the initial state in a reaction problem has amplitude x a in channel a, then 



*(all r a ->• oo) = ^A a 



{x a Fi a (r] a ,k a r a ) +y a Gi a (f] a ,k a r a )} /r a . (12) 

J 



In this expression, channel a is characterized by daughter nuclei with wave functions i/'af 1 ano - ' l l } a2 2 i way e number 
k a , Sommerfeld parameter r] a , and daughter separation r a . The orbital angular momentum of this channel is l a , 
and the square brackets denote coupling of the daughter angular momenta J a2 and l a first to "jj-coupled" angular 
momentum j a and then with J a2 to form total angular momentum J. (The use of this coupling anticipates its later 
utility in defining channels in QMC wave functions.) The antisymmetrization operator A a carries out an antisymmetric 
sum over all partitions of nucleons into daughters al and a2 with mass numbers A\ and A 2 and multiplies by the 
normalization A\\A 2 \/A\. The index a is taken to specify the daughter nuclei as well as the channel quantum 
numbers j a , J al , J a2 , l a , and ir a . 

Away from resonance, the irregular-function amplitudes are given by y a = Ylb Kab^b- At a pole of K a i,, all 
asymptotic channels that are coupled to the resonance have wave functions proportional to Gi a {rj a ,k a r a ). Then the 
corresponding x a are irrelevant and may be set to zero. If '3/ at a pole is normalized to unit probability within some 
finite volume that includes the whole interaction region, the residue of K aa is proportional to \y a \ 2 |44j,|45|. The partial 
width of the resonance in channel a is proportional to this residue, just as it was in the single-channel case discussed 
above (46l [48j . It can be shown that the imaginary part of the corresponding S'-matrix pole energy is proportional to 
the sum of these partial widths 48], so the pole residues give the partial widths, whereas the S'-matrix pole location 
gives only the total width. 

The channel radial functions defined by the summed terms in Eq. (1121) can be isolated from ^ by projecting it onto 
the channel functions 



$a(£al,£a2,r a ) = Ai$a,z>(£al , £a2 , r a ) 



(13) 
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where the channel function in a given partition p of the nucleons into clusters is 



^a,p(£ali £a2, r a ) 



^f(Cl)[^ 2 (^ 2 )^(fa) 



(14) 
(15) 



The £ a i are the internal coordinates of ip > written in terms of nucleon coordinates in the specified partition p as 
differences of the coordinates rj of nucleons within daughter i and the center of mass of daughter i. The daughter 
centers of mass v ai are related to their separation r a by 



r Q 2 - r Q i = r a 
Air al + A 2 r a2 = . 



(16) 
(17) 



Using bracket notation to indicate inner product in nucleon spin-isospin space and integration over all nucleon coor- 
dinates R = {ij}, 



A 1 \A 2 \ 



A\ 



which for an exact solution of the Hamiltonian gives the overlap function 

R a (r) ee (<I> a \S(r a -r)\y)/r 2 

r^oo {x a Fi a (r) a ,k a r) + y a Gi a (r]a, k a r)} /r . 



(18) 



(19) 
(20) 



It is this function that multiplies the daughter wave functions in Eq. (fT2"j). (Formally the integral extends over 
only 3(vl — 1) coordinates, because center-of-mass motion is irrelevant. The wave functions used in nuclear QMC 
calculations are translationally invariant, allowing the integral to be written as extending over all coordinates; the 
effect of center-of-mass motion cancels out when dividing computed quantities by the wave function normalization.) 



B. Integral relations and asymptotic normalizations 

Eq. (fT0|) suggests that the partial width of a real-energy resonance state may be found by first computing \P at 
resonance and then projecting onto the desired channel and dividing by Gi/r to obtain the partial width 



„ Kk„ 

r Q = — 

Ma 



R a (r)r 



Gi(r]a, k a r) 



(21) 



This approach is useful when $ can be computed accurately in the large-r region (as in Ref. [43]). However, it often 
happens that ^ is computed accurately in the interaction region, but computation in the asymptotic region described 
by Eq. (fT2j) is difficult or inconvenient. 

A more robust approach proceeds through a Green's function formalism. One begins with the Schrodinger equation, 

(H - E)V = 0, (22) 

with H the Hamiltonian operator and E the energy. To isolate a particular channel a, as defined above, it is useful 
to choose a partition p of nucleons into the daughter nuclei in that channel and write H as a sum of parts internal to 
the daughters and parts that depend on their relative motion. Working in the center-of-mass frame, 

H = T a , p + [T'f + H^ p + H^ p + V§- Vg, (23) 

where T a ^ p is the kinetic energy of relative motion of the daughters. The operator H®' p is the part of the Hamiltonian 
(kinetic plus potential) involving only the coordinates within daughter i. U^f contains the remaining terms of 
the potential, consisting of interactions between nucleons in daughter 1 and nucleons in daughter 2: 

E£f= E + \ E v ^ + l E Vm, (24) 

iGal;i£a2 iEal;j,kEa2 ij£al;k£a2 
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This is not an effective interaction but rather the sum of all terms of the many-body potential (with two-body 
terms Wy and three-body terms Vijf.) linking the two daughters. The point-Coulomb interaction between daughters, 
YS( r a) = ZaiZ a 2e 2 /r a , is added and subtracted for reasons that will become apparent later; it is zero if one of the 
daughters is a neutron. 

The energy E may be similarly broken up into the sum of the daughter internal energies E a i and the channel energy 
E a : 



E — E al + E a2 + E a . 

Then Eq. ([22]) becomes 

( T a, P + tCf + VS-VS + HZ* + - E aX - E a2 - E a )f = 0. 
Rearranging terms and multiplying by the operator [T a p + Vq — Ea} -1 , we find 

* = -[T a , P + Vg - E a ]- X {U:« - VS)* 

- {T a , P + VS- E^^ {Ht' P + H%' p - E al - E a2 ) * . 



(25) 
(26) 
(27) 



The next step is to project onto channel a, as in Eq. (fP9|) . Projection of Eq. (|27[) onto a channel function $ a;P 
of Eq. ([TJ} gives 



*a,p(r'Jj [T a<p + VS E a ]- L (U^ - VS) * 

$ a , p (r' Q )l f [T a , p + VS - E^- 1 (H^ p + H^ p - E al - E a2 ) * 



(28) 



(abbreviating $ a p on the right side by omitting the cluster internal coordinates), but since 



the second line of Eq. (|28|) is zero, and 



* = - 



Kp(0 ' [Ta.p + VS - Ea}' 1 (U% - VS) * 



(29) 



Now, the operator \T a ,p + V§ — E a ] 1 takes functions of nucleon coordinates R to functions of nucleon coordinates 
R' with different values of the separation r a but with the £f untouched. Its application to a function </>i(R) and 
projection onto a second function <p2 (R') may be written as an integral over a Green's function that contains a product 
of Coulomb wave functions: 



4(R') [T a , p +VS- Ea}' 1 <MR) = j^4(R') / dPr a 



3 i 7 i a (»/a,fca»'<)G'; a (?7 a ,fc Q r > ) 



r<r> 



YiMYfMM*), (so) 



following the usual notation that r< denotes the smaller of r a (a Jacobi coordinate specified for partition p by R) 
and r' a (specified by p and R'), while r> denotes the larger of r a and r' a . Rewriting Eq. (TTS]) in terms of Eq. (j2"9")l and 
applying Eq. (|30p , the result of integration over r' a and antisymmetrization over partitions p is 



<* a |^V^l*} = -> 



h?k n 



Fi a {r]a, k a r) 



•A-a 

r a <r 



Fi a {Va, k a r a ) 



$a, P (r a ) (U^-VS)^d iA R 



(31) 



Gljj]a,k a ra) 



$a, P (r Q ) (U T a *-VS)*d 3A R 



The constant of integration B a oo can be determined from 4" as in Ref. [36|, yaj; as discussed below, B a oo = at 
resonance. Finally, I denote the integrals in Eq. (|3"Tj) by 



Ba(r) 
C a (r) 



2^, 

~Vka 

Wk a 



Gi a (Va,k a r a ) 



Fi a {Va, k a r a ) 



(U^-VS)*d 3A R 
*aAra)]\u^-VS)*d 3A R, 



(32) 
(33) 
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so that 

R a {r) =-{[B a , 00 +B a (r)]Fi(r ja ,k a r)+C a (r)G l (r ]a ,k a r)}/r. (34) 

Many derivations of Eq. (I34[) and equivalent (through analytic continuation) bound-state expressions may be found 
in the literature (e.g. Refs. j36l.l47l.l49l.[50j). 

Several observations may be made at this point. The choice to add and subtract the point-Coulomb interaction Vq 
in Eq. (l23|) has two important consequences. First, it guarantees that the overlap R a {r — > oo) computed from Eq. ([34]) 
is a linear combination of Fi and Gi, as required for solutions of the Schrodinger equation governing ^ . Second, the 
nuclear interaction is short ranged, so that at large separations, the interaction between daughter nuclei is dominated 
by the monopole term of their Coulomb interaction. Since this is just equal to Vq, the difference (f/^f — Vq) goes 
rapidly to zero for r beyond the range of the nuclear interaction. For the interactions and wave functions discussed 
below, this typically happens at r ~ 7 fm. Thus, B a (r > 7 fm) ~ 0, and C a (r) takes on approximately its asymptotic 
value for all r > 7 fm. These two properties make Eq. (|34|l especially useful for computing asymptotic properties of 
ty. Indeed, similar integrals tend to appear in scattering theory for exactly this reason. 

Comparison of Eqs. (|34"]) and ^ shows that since B a (r — » oo) = 0, 

+ x is g a (r -> oo) 

tand a = A aa = , (35) 



as discussed in many places in the literature (e.g. [3J, |36j, |48( ) . In fact, Eq. ff35|) may be improved upon significantly to 
yield estimates of K aa that are second-order variational in approximations to VP [la, [H, 13, [EH . (These improvements 
require considerably more computation, so I do not pursue them here.) From the discussion in Sec. Ill Al at resonance 
there is no contribution from the regular function so that B a ^ = and K aa has a pole, 
pole (for "J normalized to unity within some radius r norm that contains the interaction region) is 



K = ^\Ca(r -> oo)| 2 , (36) 



up to corrections of order 



n | ^norin-^o (^norm) | 7 
2^a dE a 



Gi a (rj a , fc a r norm ) 



(37) 



discussed in Refs. Q HE} . It is easy to show from Eq. ((SJ that the residue of the 5-matrix pole is twice the residue 
of the corresponding if -matrix pole, and thus the partial width of the resonance in channel a is given by 

r a = — |C a (r^oo)| 2 , (38) 

Ma 

again up to corrections of order e. Since the integral defining C a (Eq. (|3"3"|) ) is short-ranged, Eq. (|3"5)) for the partial 
width may be formulated as the square of a straightforward integral over all nucleon coordinates without specifying 
a boundary. 

This result for the width is easily generalized to include bound states. These occur at negative energy, E a = —\E a \, 
so that k a = i\k a \ and rj a — —i\i] a \. Since they must be square-integrable, their overlap functions are asymptotically 

R a {r ^ oo) oc Oi{ria,k a r ->■ oo)/r oc e ^ r -^-x, (39) 

For purely imaginary k a , Oi is proportional to the Whittaker function, W_^ a ^ l+ i(2\k a \r). Working with analytic 
continuations of the linearly-independent pair of functions Fi and Oi instead of F\ and Gi , an integral relation analogous 
to Eq. ([3T]) may be derived. It yields a result for R a {r) that is guaranteed to have the correct form (Eq. (|39l) ) at large 
r even if it is generated from an approximate ^ that does not solve the Schrodinger equation exactly. In fact, an early 
use of integral relations of the form considered here was to produce overlaps with the correct asymptotic behavior 
from asymptotically- incorrect Hartree-Fock wave functions [43 . |52| . They have been used more recently to compute 
the asymptotic normalizations a a defined by 

R a (r -> oo) = a a W_ lvMa+i (2\k a \r)/r (40) 



from shell- model and ah initio wave functions [3CH33L l50l l53l . |54| . 

It is evident from the discussion above that a variational approximation to the wave function ^ allows two calcu- 
lations of R a (r): one from the definition in Eq. (|T5)) and one from the integral relation of Eq. (|3~Tj) . If "P, ip^i, and 
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ip a 2 are all eigenstates of their respective Hamiltonians, both methods yield the same result. In a typical application 
of many-body methods, \& is quite accurate in the interaction region but inaccurate where the daughter nuclei are 
widely separated. Eq. ([TBI) yields R a (r) that is only as accurate as f is at r and is not guaranteed to have the 
correct form at large r. If R a is computed from Eq. (|31[) . its value at any radius depends only on the values of W 
within the interaction region (apart from the integration constant -B a ,oo, which is nonzero only for nonresonant open 
channels). At large radius the correct asymptotic shape is guaranteed. It seems likely that since Eq. (|3Tj) incorporates 
the Hamiltonian directly (through ?7™M and not just through the approximate 'J, it is also more accurate at smaller 
r than the directly-computed R a (r) [521 ] - 

III. APPLICATION TO VARIATIONAL WAVE FUNCTIONS 

In the following calculations I use variational wave functions that were computed by the VMC method. They 
are approximate solutions for a Hamiltonian consisting of the sum of the Argonne vis (AV18) two-body [55[ and 
Urbana IX (UIX) three-body [5(| interactions, which appear both in the computation of the wave functions and in 
the operator U^ p used in the integral relation. In this section I describe enough of the structure of the VMC p-shell 
wave functions to discuss their asymptotic properties, I explain the implementation of the integral relation, and I 
finally comment on the asymptotics of the wave functions. The wave functions are described in detail in Ref. [57| . 

A. Variational wave functions 

The VMC wave functions begin with the spin-isospin function $ Q (0000) 1234, in which the spins and isospins of the 
first four particles are organized into a Slater determinant like those of a filled, a-particle-like 0s shell. (The zeroes 
denote the total angular momentum </, its projection M, the total isospin T, and its projection T z ; this notation 
follows conventions in previous QMC work and should not be confused with the channel function $ a ). The addition of 
further particles while retaining antisymmetry requires spatial dependence in the wave function, i.e., placing particles 
into the p-shell. This is done using orbitals </>p S (r a i) defined in terms of the separation r a i of particle i from the 
center of mass of the "alpha core." The orbitals are derived from Woods-Saxon potential wells and coupled to form 
states of definite angular momentum, parity, spatial symmetry, and isospin in an LS'-coupled basis, and they fall off 
exponentially at large r a i. The effects of pairwise interactions between nucleons (through the AV18 potential) are then 
accounted for using two-body scalar correlations f(rij) that mainly account for particles' avoidance of the repulsive 
core of the potential and are functions of the separation nj of particles i and j. It has been found advantageous to 
allow different scalar correlations depending on whether particles i and j are both in the s-shell (f ss ), both in the 
p-shell (f pp ), or one in each (f sp ). Finally, there are analogous three-body scalar correlations (f sss , f spp ) due to 
both two- and three-body interactions. All of these elements (a core, LS coupled orbitals, and scalar correlations) 
are antisymmetrized among particle permutations to make up the Jastrow wave function: 

*j =m n fm n n n n f ss ^ 

\ i<j<k<4 n<4 5<m<A m<q<A t<u<4 

xj] n / sp m n / pp (^) 

i<4 5<j<A 5<k<l<A 

x (0LS[n]®A(LS[n]JMTT z )p) I . (41) 

LS[n] J 

The index P denotes a specific permutation of the particles into s- and p-shells (subsequently antisymmetrized by the 
operator A) and the amplitudes PlsIu] specify the admixtures of p-shell LS states of Young diagram [n]. In turn, the 
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configuration of good L, S, [n], J, M, T, and T z is given by 
$ A (LSJM[n]TT z ) P = $ a (00 00) 1234 J] 



5<i<A 



5<j<A 



LM L 



Jl Xk{\m k ) 



5<k<A 



SM S 



JM 



n «4*« 



5<1<A 



(42) 



TT, 



where the spinors x% an d v % specify the angular momentum and isospin states of particle i. Because the Hamiltonian 
includes operator as well as scalar terms (i.e., it acts on particle spins and isospins), there are operator correlations 
in addition to central correlations. These are accounted for by writing the variational trial function as 



S\{\l + U K 



(43) 



where Uij and U™ 1 are two- and three-body operators, and S is a symmetrization operator, needed to preserve the 

ry o 

of the form 



antisymmetry of ^ j because the Uij and U™ 1 do not commute amongst themselves. The operator correlations are 



2<<}<6 



U q(nj)0?j 



(44) 



where is an operator-independent three-body correlation, and the operators 0\j — T{ ■ Tj, <Ji ■ tjj, <Ji ■ (TjTi ■ Tj, 
Sij, and SijTi ■ Tj (where <Ji are nucleon spin operators, Tj, are isospin operators, and Sij is the tensor operator), 
appear in the largest operator terms in the AV18 potential. The w 9 (r^), together with the scalar correlation f ss (rij), 
solve a set of coupled Euler-Lagrange equations with coefficients that serve as variational parameters, discussed in 
Rcf. [Hj]. The central correlation f ss falls off exponentially to reflect the strong binding of the s-shell particles, as do 
the u q functions. 

The f sp (rij) and f pp {rij) correlations are constructed to approximate f ss (rij) at small but to approach constant 
values at large Tij . This guarantees that where particles approach each other closely the wave function is governed by 
the nucleon-nucleon interaction, but that the correlation between widely-separated particles is dominated by binding 
to a "mean field" accounted for in the 4>p S ^ orbitals. Thus, the asymptotic region of is dominated by the </>p S '™', 
which have much longer tails than the f ss . 

All of the functions appearing so far in this section are specified as functions of variational parameters, either 
explicitly or in the differential equations solved to compute correlations. The optimum values of those parameters are 
found by searching the parameter space to minimize the energy expectation value, 



E < 



(*t|*t) 



(45) 



with both numerator and denominator computed by Monte Carlo integrals over the particle coordinates. As a final 
step, the Hamiltonian is diagonalized with respect to the p-shell configurations labeled by LS[n] to find the coefficients 
PlsM- Diagonalization both improves the variational energy of the ground state of given J T and T and provides 
access to excited states. 

Energies of unbound p-shell states can almost always be lowered by making their wave functions more diffuse 
(closer to threshold). This is also often true of bound states, where variational energies (but not GFMC energies) 
can lie artificially above breakup thresholds because of shortcomings of the variational ansatz. For both bound and 
resonance states, this problem is addressed by constraining variational parameter searches to keep charge radii close 
to experimentally known ground-state charge radii. 

The values of the variational parameters for all states used in the present calculations were provided by R. B. Wiringa 
[59] | . They are the results of calculations in bases of good isospin in which individual nucleons typically cannot have 
definite charges. For convenience in defining the integral relations for neutron or proton decays, I carry out the 
calculations below in a basis of good nucleon charge (essentially an m-scheme in the particle isospins) so that the 
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emitted nucleon is definitely either a neutron or a proton. For given variational parameters, this is only a change of 
representation and does not alter observables. The variational minimization has in all cases been carried out for the 
state of lowest T z in each isomultiplet. I compute widths of the isobaric analogues of these states by using isospin 
rotations of the minimized states rather than carrying out independent variational minimizations. 



B. Overlap and integral-relation calculations 

Calculations of explicit overlap functions R a (r) (Eq. (jT9j) ) and integral relations (Eqs. (|32|) and (J33J) ) involve 
several common elements. While integral relations are much more time-consuming because they contain the opera- 
tor £7 r "f , most of the computational tasks in organizing the calculation amount to constructing the channel vector 
&a,p(£ai> £a2> r a) gi yen by Eq. (fT5|) and then contracting either or (U^[ — Vc)^t against it. For nucleon emission, 
I take "^ai 1 to be the wave function of the daughter nucleus and - 0a2 2 to be the spinor of the emitted nucleon. 

The integrals of Eqs. (fT8|) and (|33|) are computed by Monte Carlo integration over the particle coordinates, using 
the same sampling algorithm that has long been used to compute the energy expectation value, Eq. (|45j). Sampling 
follows the Metropolis algorithm, using the weight function 

W(R) = *£(R)# T (B.) , (46) 

where \ty(R) is the variational wave function for the A-body parent nucleus at particle coordinates R. The delta 
function in Eq. (|18p is accounted for by sampling all particle coordinates and sorting the samples into narrow bins of 
specified channel radius; this builds up the entire function R a (r) from a single Monte Carlo walk. The normalization 
integral needed to give ipai 1 Uln t n orm is computed in the same Monte Carlo walk by which the overlap or integral 
relation is computed. 

Only relative coordinates are used in the definitions of and $ a ,p, so the results of the calculations are all 
explicitly translation invariant, and for operator 

. . _ J S(r a — r)/r^ spectroscopic overlap 



the quantity computed is 



Fi(r] a ,k a r a ){U^ - V§)/r a integral relation 
(^ a \M\^ T ) 



(48) 



The routines used to compute the integral relations were written as modified versions of existing spectroscopic-overlap 
routines [60l - l63l |. The integral- relation routines were used previously to compute bound-state ANCs 33], and only 
very minor modification (replacing regular Whittaker functions with regular scattering functions) was necessary for 
width calculations. 

Finally, the operator U^f must be considered. It is just the potential-energy operator of the AV18+UIX Hamil- 
tonian, but with the restriction that only terms involving the p th (emitted) nucleon are considered so that Eq. (|24|) 
becomes: 

tcr = E^+ E v *ip- ( 49 ) 

Its action on is evaluated by calling the potential-energy routines with instructions to omit all terms purely internal 
to ip^i 1 f° r * ne gi yen permutation. It is then straightforward to contract (C/^'f ~ Vc)^t with [$o,p(£oX> £a2> r a)V f° r 
a given configuration R. 

It is instructive to examine the integrand of Eq. (|33j) by inserting a delta function 5{r a — r)/r^ into the integral. 
This is evaluated just like the delta function in Eq. (|T5| . by carrying out the Monte Carlo walk for the full integral and 
binning the Monte Carlo samples according to the channel radius instead of summing them all together. A typical 
result is shown in Fig. [lj and nearly all cases that I have computed appear similar to this graph apart from an overall 
scaling. It is evident that the largest contribution comes from r a ps 2.5 fm. It is also evident that the form of the 
VMC wave function beyond ~ 7 fm is irrelevant for the width, because the factor (U^\ — Vg) is zero there. This 
limits the integral to smaller radii and guarantees convergence of the integral regardless of what lies in the tails of the 
variational wave function. 

The motivation presented above for the integral relations assumes that ^t, "^ai 1 1 an< ^ ^ia are a ^ energy eigenstates. 
To the extent that the VMC wave functions approximate these states, the integral relations approximate the overlaps 
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r (fm) 

FIG. 1: (Color online) The distribution of the integrand of Eq. (|33[) in two j a channels is shown as a function of channel radius 
r after integrating over all other coordinates. The case shown is the 9 Li(3/2^") overlap onto 8 Li(g.s.) + neutron, but the shapes 
of the curves are similar in nearly all calculations from VMC wave functions. Discrete points show the integrand dC a /dr. The 
solid curve shows its cumulative integral, C a (r), which at large radius is equal to the wave function asymptotic normalization. 
The integrand goes to zero at small radius because it contains the regular function Fi(rj a ,k a r a ) and at r a ~ 7 fm because 
UrA —Vq—^O there. The light dotted curve indicates the distribution of samples in the Monte Carlo integration (with zero 
at the bottom of the graph). 

and widths of the Hamiltonian for which they are approximate solutions. Two difficulties then present themselves 
in applying the integral relations to VMC wave functions: 1) comparison of the results to other calculations using 
the AV18+UIX Hamiltonian is not possible because (except for 5 He) no such calculations have been done by other 
methods, and 2) comparison to experiment is complicated because experimental resonance energies are not reproduced 
exactly by the Hamiltonian. 

The mismatch of resonance energies from experiment can be dealt with straightforwardly. The integral relations 
of Eqs. (|32p and (|33[) require an assumed channel energy E a , from which k a and r\ a are computed. Formally E a 
should be the channel energy of the AV18+UIX Hamiltonian (known from GFMC), but the mismatch between this 
and the experimental E a is often ~ 1 McV. Since the potential energy in the p-shell is much larger (—100 to —400 
MeV), the experimental E a is close enough to the AV18+UIX channel energy that it can plausibly be used in the 
integrals. (See Sec. II V AD This choice accounts for most of the well-known strong dependence of the width on the 
resonance energy so that comparison with experiment is possible. It was found previously that using the experimental 
separation energy in integral-relation calculations of bound-state ANCs produces results in generally good agreement 
with experiment [H| . Similarly, the results below indicate that using the experimentally-measured resonance energy 
as E a allows prediction of experimental widths. I present the results of using both experimental and GFMC values 
of E a in the integral relations. 

More consistent calculations will require Hamiltonians that more precisely reproduce thresholds and resonance 
energies. Such Hamiltonians exist in the form of the Illinois three-body potentials pi, HH , but they have not yet been 
incorporated into the VMC code used here, and U^f for these interactions requires considerably more computation 
than for the UIX interaction. 



C. The asymptotic forms of VMC wave functions 

Consider configurations in which particle A (before antisymmetrization) is far from the first A— 1 particles. Because 
f pp and f sp approach constants and f 3S decays rapidly at large r^, the shape of *Bt in this part of the wave function 



12 



is dominated by the shapes of the single-particle functions (j) p ' n| . This might be expected to give 

R a (r ^oo)^Yl lLS[n]4>n S[n] ("LS[n]r) (50) 
LS[n] 

for some amplitudes "fzsin] and Jacobian-like factors u)Ls[ n ] that emerge from the correlations of Eqs. (|4lT) -(j43 |) when 
the overlap integral is computed. The factors u>LS[n] account for the distinction between the channel radius r and 
the distance r a i of a p-shell particle from the center of mass of the s-shell core; these only coincide when A = 5 and 
otherwise differ by the mean difference between nucleon-core and nucleon-daughter distances. 

In general, the cj)p S ^ that emerge from the variational procedure do not yield the correct long-range asymptotic 
shapes for the overlaps R a . This is most readily seen for open channels, where solutions that should oscillate at 
large channel radius instead fall off with an assumed exponential dependence. In closed channels, the condition of 
square integrability gives zero for the analogue of B atOC in Eqs. (j3"Tj) and so that Eq. (HOI) holds for the true 

wave function. Because no single LS[n] term typically dominates a given i? a (r), it is in general difficult to construct 
^LS[n] ^ o ga y g fy pU|) for all possible channels at once. This difficulty is compounded by the problem that the 
energy expectation value driving the variational minimization receives very little contribution from the wave function 
tails, so the variational principle does not constrain the low-probability tails of the wave function very strongly. The 
application of the integral approach to bound state ANCs in Ref. [33| avoided these difficulties and effectively matched 
the correct asymptotic form onto the better-computed interior of the wave function. 

Pseudobound VMC wave functions approximate resonance states with square-integrable wave functions, in which 
the f SB and functions cut off the wave function tails exponentially. This exponential cutoff can be understood 

at "medium" range (4-8 fm) as forcing a resonance form on the wave function and beyond this range as providing a 
regularization to normalize the unbound state despite its formally nonzero amplitude at large radius. Such a regu- 
larization is important for quantities like electromagnetic transition strengths (e.g. (65l . I66l0 and for the approximate 
relation between the asymptotic normalization and the partial width given in Eq. (|38p [45j . As long as the resonance 
wave function is computed reasonably accurately within the region where the integral of Eq. (|33p is nonzero, and it 
is normalized to unity over the region where \E't is larger in amplitude than in the asymptotic region, the integral 
approximates the asymptotic normalization of Gi and thus the partial width. The cutoff of at large nucleon 
separation is illustrated with overlap functions in Fig. [2 

IV. RESULTS 

I now apply the integral relation of Eqs. p8[) and Q33p to compute widths of several unbound energy levels in 
nuclei of mass numbers 5 < A < 9, from VMC wave functions. I mostly choose energy levels that empirically have 
small width (under 1 MeV) and are dominated by nucleon emission. The integral relation is valid for decay channels 
in which all three wave functions (one parent and two daughters) are known, but I limit this first examination to 
nucleon-emission channels, where one "daughter nucleus" is a neutron or proton. I concentrate on two-body final 
states, but I also model some three-body decays as sequential processes, e.g. the 9 B ground state decaying to a proton 
and the unbound ground state of 8 Be. The subsequent decay of 8 Be to two alpha particles can be neglected because 
of its small width. (In principle the widths of unbound daughter states should be integrated over, but I find that in 
all cases considered this correction is much smaller than either the experimental errors or the widths of omitted decay 
channels.) 

The results of the width calculations are shown in Table HI Each calculation was carried out twice: once assuming 
the channel energy E a from the AV18+UIX Hamiltonian (known from GFMC calculations) and once using the 
experimentally-known resonance energy. (For the second J n = 2 + states in 8 B and 8 Li and the J* = 7/2~ states in 
9 Be and 9 B, AV18+UIX energies have never been computed with GFMC). Where the experimental channel energy 
is unknown or uncertain, I have used instead the GFMC energy for the AV18+IL7 Hamiltonian, which gives a better 
overall fit to experimental energies than AV18+UIX. 

I have taken the experimental energies and widths in Table U mainly from data compilations [67], [68| . Nearly all of 
the widths in the compilations are "observed" widths, which coincide with the full width at half maximum (FWHM) 
of cross section peaks and are proportional to the sum of the if-matrix pole residues. The integral-relation widths 
should correspond to these quantities. In two cases ( 7 Li(5/2^~) and 9 Be(l/2 - )), examination of the source literature 
indicated that the available experimental widths are "formal" widths of the i?-matrix formalism. I have converted 
the experimental widths listed for those states to "observed" widths by the usual relation 



(51) 
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t — — i — r| — i — i — i — | — i — i — i — | — i — i — i — | — i — i — r 




r (fm) 

FIG. 2: (Color online) Overlaps R a (r) of the unbound J n — 1 + state of 8 B onto 7 Be + proton, computed directly from Eq. (|18p 
(discrete points with Monte Carlo error bars) and from the integral relation of Eq. pip (solid curves). At large radius, the 
directly-computed overlap falls off exponentially because the (f>p Sln] in the VMC wave function fall off exponentially. However, 
the integral-relation overlap (with B a .ao = 0) gives the exact proportionality to Gi expected at resonance. The absolute 
normalization Cf^ = C a (r — ¥ oo) multiplying Gi gives the partial width of 8 B(1 + ) in each 7 Be + p channel; the normalized 
asymptotic is shown as a dashed curve for each channel and manifestly merges with the full overlap at > 4 fm. 



where b (taken to be 4 fm) is a matching radius, Si(E) is the shift function of i?-matrix theory [40|, prime denotes its 
derivative, and 7^ is defined by the formal width T a = 2j%Pi„(k a b) and the penetrability function Pi a (k a r). (See the 
section "Definitions of resonance parameters" in Ref. [67|.) The sources of the experimental widths are indicated in 
Table i 

The widths presented in the table are sums over all p- and /-wave decay channels. The predicted /-wave contribu- 
tions are less than 1% of the total in all cases except the decay of 9 Li(7/2~) to 8 Li(2 + ), where it is computed to be 
23% of the total. The p-wetve decays are in most cases sums of P1/2 and P3/2 contributions. The decompositions into 
P1/2 and P3/2 (or into a channel-spin coupling scheme) are available on request to the author. 

The table only includes cases for which established experimental or GFMC energies are available. I have also 
computed widths of several states of 9 He with varying assumptions about resonance energies as described below, but 
I have not included these numbers in Table [IJ 

The computed widths of those states and channels for which experimental data are available are shown graphically in 
Fig. [31 Only for those states where partial widths are available from experiment do I show partial widths in the graphs 
- otherwise total widths are shown. Fig. [3] demonstrates the wide dynamic range of the integral method, extending 
from 0.0005 to roughly 1 MeV. With the exception only of the very broad 5 He states, which present problems 
discussed below, all computed widths are within a factor of three of experiment - within a factor of two if states 
with uncomputed alpha and direct three-body width are omitted. The error-weighted mean ratio (r; n tegrai/r e xpt) of 
the integral-method width to the experimental width is 0.82 ± 0.29, and the x 2 statistic for the difference between 
computed and experimental widths is 5.9 per degree of freedom (omitting 5 He for a total of 19 states). Restricting 
consideration to those states with no omitted channels gives (rj n t egra i/r oX pt) = 1.09 ± 0.04 and xt — 1-5 with v = 9. 
The errors used to compute these statistics are mainly in the measured width r cxp t. 

Setting aside the resonance energies, the main sources of uncertainty in the theoretical calculation are in the 
accuracy of the variational wave function and of the potential. Because it is unclear how to estimate these errors, the 
theoretical errors reported in Table U are the quadrature sums of Monte Carlo sampling errors with errors propagated 
from the input resonance energies. Errors propagated from experimental energies are typically small compared with 
the other errors. However, each GFMC channel energy is the difference of a resonance energy and a threshold energy 
from separate calculations. The error on this difference is often large compared with the channel energy, and this 
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TABLE I: The results of integral-relation calculations of widths. Results are shown from calculations in which the channel 
energy was assumed equal to its experimental value ("From Expt energy") and to the AV18+UIX channel energy known from 
GFMC ("From AV18+UIX energy"). Error estimates are quadrature sums of Monte Carlo sampling errors and uncertainties 
from the input energies. Where no experimental energy is available, the results in the "Experiment" columns were computed 
using the GFMC energy with the AV18+IL7 Hamiltonian, and they are indicated by square brackets. The column "Matches 
7r/2?" indicates whether the overlap function appears consistent with a resonance state, as discussed in Sec. II V Al See Sec. IIV Al 
and Eq. (|52p for the definition of £. Energies are relative to the decay threshold in the center-of-mass frame. Experimental 
energies are taken from data compilations 67, 68] unless otherwise noted. 



State 



Daughter Experiment 
E (MeV) T (MeV) 



From Expt energy From AV18+UIX energy Matches £ 
Tvmc (MeV) £uix (MeV) Tvmc (MeV) tt/2? 



°He(3/2-) 4 He(0 + ) 
5 He(l/2-) 4 He(0+) 
7 He(3/2") 6 He(0+) 


0.798 
2.07 

0.445(3) 


0.648 [67] 
5.57 [67] 
0.122(13)" 


0.307(5) 

0.582(13) 

0.114(12) 


1.39 
2.4 

1.68(13) 


0.684(11) 
0.711(15) 
0.77(10) 


no 
no 

yes 


0.460 
0.429 
0.092 


7 He(l/2") 6 He(0+) 
7 He(l/2-) 6 He(2+) 
7 He(l/2") sum 


3.05(10)" 
1.25(10)" 


2.0(1.0) c 


1.98(12) 

0.42(6) 

2.40(12) d 


2.83(13) 
0.89(13) 
2.83(13) 


1.80(12) 

0.26(5) 

2.22(ll) d 


no 

yes 


0.21 
0.067 


Y He(5/2") b He(2+) 
7 Li(5/2 2 -) 6 Li(l+) 
7 Be(5/2 2 -) 6 Li(l+) 


1.55(3)" 

0.204 

1.60(6) 


1.99(17) [67] 
0.065(3) [67] 
0.19(5) [69] 


1.29(12) d 

0.0483(17) d 

0.43(4) d 


1.87(13) 
1.57(17) 
2.65(17) 


1.7(2) d 

0.92(13) d 

l.ll(14) d 


no 

yes 
yes 


0.165 
0.055 
0.055 


s Li(3+) Y Li(3/2~) 
8 Li(0+) 7 Li(3/2~) 
8 Li(0+) 7 Li(l/2~) 


0.223(3) 

[0.97(13)] 

[0.847(14)] 


0.032(3) [70] 


0.0344(18) 

[0.37(7)] 

[0.81(2)] 


2.10(18) 
0.67(17) 
0.78(17) 


0.88(11) 

0.24(8) 

0.7(2) 


yes 
no 
no 


0.007 
0.005 
0.004 


B Li(2+) y Li(3/2-) 
8 Li(2+) 7 Li(l/2~) 
8 Be(l+) T=V 7 Li(3/2-) 
8 Be(l+) T = e 7 Li(3/2") 
8 Be(l+) sum 6 7 Li(3/2") 


[2.18(16)] 
[2.06(19)] 
0.385(1) 
0.895(5) 


0.149(6) [68] 


[l.oo(ii)] 

[0.105(14)] 
0.0089(3) 
0.152(4) 
0.161(4) 


1.26(19) 
0.51(21) 


0.17(5) 
04+ 05 
0.21(5) 


yes 
yes 
yes 
yes 
yes 


0.004 
0.010 
0.003 
0.003 


8 Be(3 + ) T = V Y Li(3/2-) 
8 Be(3+) T — e 7 Li(3/2~) 
8 Be(3+) T = V 7 Be(3/2"; 
8 Be(3+) T = e 7 Be(3/2"; 
8 Be(3 + ) sum 6 sum 


1.81(3) 
1.98(1) 
0.17(3) 
0.335(10) 


0.50(3) [68] 


0.166(9) 

0.314(14) 

0.012(3) 

0.050(3) 

0.542(17) 


3.68(18) 
2.3(2) 
2.14(18) 
0.8(2) 


0.60(7) 
0.43(7) 
0.45(5) 
0.18(6) 
1.66(13) 


yes 
yes 
yes 
yes 
yes 


0.007 
0.003 
0.007 
0.004 


8 B(1+) Y Be(3/2-; 
8 B(1+) 7 Be(l/2-; 
8 B(1+) sum 


0.632(3) 
0.203(3) 


0.0357(6) [68] 


0.0382(15) 
0.00105(8) 
0.0394(15) 


1.3(3) 
1.4(2) 


0.26(13) 

0.5(2) 

0.8(2) 


yes 
yes 
yes 


0.001 
0.003 



8 B(3+) 
8 B(0+) 
8 B(0+) 
8 B(2+) 
8 B(2+) 
9 Li(5/2" 



7 Be(3/2 
7 Be(3/2" 
7 Be(l/2~ 
7 Be(3/2- 
7 Be(l/2- 
8 Li(2+) 



2.18(2) 
[2.55(13)] 
[2.44(14)] 
2.41(2) [71] 
1.98(2) [71] 
0.232(15) 



0.39(4/ 



0.12(4) [71] 
0.24(11) [71] 
0.10(3) [63J 



0.38(2) d 
[0.65(7)] 
[1.46(18)] 
0.51(2) 
0.039(2) 
0.145(14) 



3.7(2) 

2.17(17) 

2.30(13) 



1.08(11) 

0.47(8) 

1.3(2) 



0.98(44) 1.2(8) 



9 B(7/2" 
9 C(l/2~ 



8 Be(2+) 
8 B(2+) 



4.13(6) 
0.918(11) 



2.0(2) [68] 
0.10(2) [76J 



0.82(4) d 
0.102(5) 



1.5(3) 



0.43(26) 



yes 
no 
no 
yes 
yes 
yes 



yes 
yes 



0.007 
0.005 
0.004 
0.004 
0.010 
0.003 



a Li(7/2") 9 
9 Li(7/2") 9 
9 Li(7/2-)» 


8 Li(2+) 
8 Li(3+) 
sum 


2.366(15) 
0.111(15) 


0.04(2) [63] 


0.0012(5) 

0.043(8) 

0.044(8) d 


3.6(3) 
0.23(35) 


0.0031(13) 

< 0.50 

< 0.50 d 


no 

yes 


0.045 
0.006 


9 Li(3/2 2 -)« 
9 Li(3/2 2 )° 
9 Li(3/2 2 )« 


8 Li(2+) 
8 Li(l+) 
sum 


1.32(6) 
0.34(6) 


0.6(1) [63] 


0.52(4) 
0.17(5) 
0.69(6) d 


1.5(4) 
0.5(4) 


0.6(2) 
< 0.7 
0.9(4) d 


no 

yes 


0.014 
0.006 


9 Be(l/2-) 


s Be(0+) 


1.11(12) 


0.86(g)' 1 f72] 


0.80(12) d 


4.4(6) 


4.8(8) d 


yes 


0.0005 


a Be(7/2") 
9 Be(7/2-) 
9 Be(7/2-) 


8 Be(0+) 
8 Be(2+) 
sum 


4.72(6) 
1.69(6) 


1.2(2) [73] 


0.0082(3) 

0.40(3) 

0.41(3) d 






yes 
yes 
yes 


0.005 
0.003 



0.003 
0.006 



"I have computed an "observed" width of 112 ± 15 keV from the i?-matrix formal width of Ref. [7711 and averaged it with the 160 zt 30 
keV FWHM of Ref. @. 
"From 79], based on ground state energy from |67|| . 

c This is reported in Ref. [62| as "« 2 MeV" with no quantitative error; a 1 MeV error is used in Figs. [3] [8] and|9]and in quoted statistics. 
d Open channels other than one-nucleon emission were neglected in the calculation (alpha or non-sequential). 
e See Sec. IIV SI for discussion of the effects of isospin mixing in the observed 1+ and 3 + states of 8 Be. 

^Originally reported in Ref. [80ll , this number is construed in later work as the FWHM in the lab frame; since the center-of-mass excitation 
energy is reported in the same sentence, it appears to me to be a center-of-mass width. Its error has apparently been mistranscribed in 
later references, independent of this ambiguity. 

9 Spin-parity assignments for 9 Li follow Ref. [§jj; see Sec. IIV C 31 below. 

h The data compilations [fall average the reported H-matrix formal width corresponding to this number with the much less certain 
Breit-Wigner width of Ref. [8 1| . 
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FIG. 3: (Color online) The results of Table U for those cases included in the averages. The vertical axis shows the ratio 
rintegrai/r exp t of the integral-relation width to the measured width, and the horizontal axis shows the measured width r exp t. 
The horizontal dotted line shows Integral = r Gxpt as a guide to the eye. The errors shown are dominated by those on the 
experimental widths. Experimental errors have been added in quadrature to smaller Monte Carlo statistical errors and errors 
propagated from the energy used in computing Integral- States for which important alpha or 3-body decays have been neglected 
are indicated with asterisks next to the state labels, and large (red) circles indicate overlaps that seem inconsistent with n/2 
phase shift. The outlying 5/2^" states of 7 Be and 7 Li are complicated by the need to include the small alpha decay channel in 
multichannel i?-matrix fits when extracting widths from data; this introduces a considerable spread in reported nucleon widths 



propagates to a large error on the predicted width. 



A. Correspondence of computed states to resonance states 



The results presented in Table U and Fig. Ovary in the degree to which they match the experimental widths. Some, 
like the first 3 + state of 8 Li, are close matches. Others, like the two low-lying states of 5 He, are very far from agreement 
with experiment. The 5 He cases are particularly interesting, because they are the onl y o nes for which the correct 
widths for the AV18+UIX Hamiltonian are known from explicit scattering calculations [ll|. The known AV18+UIX 
widths are 1.5 MeV for the 3/2~ state and 5.0 MeV for the l/2~ state, while applying the integral relation to the 
pseudobound VMC wave functions gives 0.68 and 0.71 MeV, respectively. The l/2~ state lies above the centrifugal 
barrier and is so broad that its phase shift does not pass through tt/2, so it is unsurprising that the integral method 
fails for this case. The 3/2~ state does not present these difficulties, so its difference from both experiment and theory 
is clearly not a shortcoming of the potential but rather of the computational methods. 

There are several reasons that the application of the integral relation to a pseudobound variational wave function 
could fail to yield the correct width for the potential. I begin by noting that for the states examined, the function 
Fi a (r] a ,k a r a ) typically does not deviate far from its leading-order dependence on k a r a over the range of r a where 
U^ei — V§ differs significantly from zero. The main effect of changing the assumed resonance energy is therefore 
to change the overall scale of the integrand in Eq. (|33j) without changing its shape much. Thus, simply using 
the experimental channel energy in the integral relation should correct rather accurately for the mismatch between 
experimental and theoretical channel energies without introducing significant distortions in the integrand. 

It is possible that the variational minimization with constrained charge radius fails to produce good approximations 
to some resonance states. It could happen that the variational ansatz is a poor match to a particular resonant 
wave function or that the VMC wave function, being only approximate, contains "contamination" from nonresonant 
continuum states. Contamination more energetic than the desired state is precisely what the GFMC method is 
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intended to remove, and much of the success of the QMC methods lies in the exclusion of low-energy excitations from 
the variational wave function. Any off-resonance contamination produces contributions to the integral relation that 
do not correspond to the pole residue. The danger of contamination would seem to grow with the resonance width, 
because a less-peaked density of states implies a wider range of states that are similar to the resonance and difficult 
to eliminate by variational minimization. The presence of nearby resonances in the same channel may compound this 
problem by contributing contamination with rather different structure from the resonant wave function being sought. 

Yet another kind of difficulty lies in the normalization of the wave function. The relationship between the width 
and the asymptotic normalization of standing-wave states found in Refs. [ID, 0, IH| depends on the wave function 
being normalized in a finite volume and its amplitude being small at the boundary of that finite volume as expressed 
in Eq. (|57|) . I have assumed that the normalization volume is effectively defined by the exponential fall-off of the f ss 

and cj)p S ^ correlations in the tails of the variational wave function. If either the true or the computed wave function 
fails to fall off rapidly enough, the normalization is problematic in ways that can be viewed either as the lack of an 
effective cutoff or as neglect of a large surface-amplitude correction. 

The normalization problem is perhaps the most straightforward to examine. The final column of Table U displays 
a parameter 



where r max is the location of the maximum |i? Q (r)|, and Ra{r) is computed directly from Eq. (|18|) . This ratio 
measures the amplitude of the variational wave function just outside the interaction region relative to that in the 
interior (accounting for the r 2 dependence of the volume of a spherical shell). We may expect to encounter difficulties 
when C > 0.1, and this occurs for four states: the two 5 He states and the l/2~ and 5/2~ states of 7 He. The relevance 
of C may also be viewed in light of Eq. (|37p . If I adopt 8 fm as the effective boundary of the normalization volume 
and assume (cf. Fig. [5] below) that typically R a (8 fm) ~ 0.02 fm~ 3 ^ 2 , then Eq. (l37l) gives e < 0.2 for all but a very 
few of the widths computed here. Much larger values of e occur for 5 He(l/2~), 7 He(l/2 _ ), and minor decay channels 
of a few other states. 

The overlap functions may also be tested for consistency with expectations for a resonance. The integral relation 
for the width is the r — > oo limit of R a as computed from Eq. (|31[) . Given a wave function, it is possible to compute 
the overlap function at all radii using Eq. (|3ip given B aj00l and the result is likely more accurate than that from direct 
calculation of R a . 

Determining B a oa from a variational wave function is a tractable problem [36l [38| . but I do not pursue it here. 
Since I have tacitly assumed in computing widths that B aao = as required for a resonance state (cf. Eq. (|35|) ). I 
can compute the full overlap function from the integral relation (Eq. (|3"Tj) ) with B atQO = and check that it matches 
the overlap function computed directly from (Eq. (|18jl ). If the two overlaps are in rough agreement in the interaction 
region, then the VMC wave function is consistent with a resonance wave function and may validly be used to compute 
a width. 

As a test of the approach, I apply it to bound states, since their integration constant corresponding to B a ao is 
zero by definition. Results are shown for several states in Fig. HJ and I have computed them for all of the channels 
considered in Ref. 33] . The agreement between the two calculations, especially for s-shcll nuclei where the VMC 
method is more accurate, is excellent. The deviations of the integral-relation overlap from the direct overlap are likely 
to be improvements: the integral relation contains more information about the potential than does the VMC wave 
function by itself, and it guarantees the correct r — > oo asymptotics. For some nuclei with A — 3, 4, and 7, GFMC 
calculations of overlaps exist (albeit for the AV18+IL7 Hamiltonian, not AV18+UIX) The results of the GFMC 
calculations (dashed curves of Fig. Q| are not severely different from those of applying the integral method to VMC 
wave functions, and they deviate from the VMC overlaps by similar amounts. This experience with bound states 
indicates that overlaps computed from integral relations are at least as accurate as directly-computed overlaps of 
VMC wave functions and are not in conflict with GFMC results. They may therefore be very useful for calculations 
of spectroscopic factors and transfer and knockout cross sections. 

Overlap functions are shown for 8 B(1 + ) in Fig. [2] and for several other resonance states in Fig. [5J These were 
computed both directly and by integral relations with B a oo — 0. In all cases, the integral relation replaces the 
artificial exponential fall-off of the VMC wave function with the oscillatory behavior of G; . Based on experience 
with bound states, the two overlaps should agree at r < 4 fm. In some cases (e.g. 9 C(l/2~)) the match there is 
quite good, while in others (e.g. (n 6 He| 7 He(l/2~))) overlaps from the two methods seem to have little to do with 
each other. I conclude that when results of the two methods are qualitatively very different inside 4 fm, there is an 
inconsistency with the assumption of ir/2 phase shift so that the pseudobound VMC wave function may not allow 
accurate width calculations. The penultimate column of Table U indicates for each width a qualitative judgment of 
whether the two methods of computing overlaps agree, and cases where they do not are indicated in Fig. [3] by large 
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FIG. 4: (Color online) Overlap functions for several bound states, computed from the bound-state analogue of the integral 
relation in Eq. (|31l) (solid curves) and from the definition in Eq. (|18[) (squares with Monte Carlo error bars.) The l a and j a 
quantum numbers of the "virtually emitted" nucleon are indicated by labels near the appropriate curves (and distinguished by 
color). Where they exist, overlaps computed directly from GFMC wave functions [60( are shown as (black) dashed curves. 

(red) circles. Examination of the table reveals that with the exception of the + state in 8 B, failures of the "7r/2 
assumption" occur only in a few neutron emission channels and in no proton emission channels. Considering only 
states consistent with 5i = tt/2, I find (ri n t e grai/r cX pt) = 0.83 ± 0.30 and xt = 6.3, this time for v = 15 instead of 
18. Further restriction to states in which all decay channels are computed gives (ri n tcgrai/r cX pt) = 1.08 ± 0.04 and 
xl — 1-6 for 8 degrees of freedom. These are essentially the same results as for the entire data set. It appears that 
the best predictor of whether a calculated width will match experiment is simply whether the calculation includes all 
channels contributing to the measured width. 

B. Isospin mixing in 8 Be 

A difficulty presents itself in considering the pair of J 71 " = 3 + states of 8 Be at ~ 19 MeV and the pair of 1 + states 
at ~ 18 MeV. Each of these doublets consists of one T = and one T = 1 state that have mixed, with the lower state 
of each pair predominantly T — 1. The VMC wave functions have definite isospin, so incorporating the mixing into 
calculations of their widths presents ambiguities. Simply ignoring the mixing and assigning the T = 1 widths to the 
lower states and T = widths to the higher states does not work. For example, in the 3 + doublet, the lower state 
has measured width 270 ± 20 keV and the upper has 227 ± 20 keV; the unmixed widths from the integral method are 
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FIG. 5: (Color online) Overlap functions of several resonance states, computed directly from the definition (Eq. ()18|l ; points 
with Monte Carlo errors) and from the integral relation of Eq. (131[) with B a ,oo = (or 8i a — n/2; solid curves). Where more 
than one angular momentum channel is available, the channels are labeled. The four cases on the left show good agreement 
between the two methods at small separation (r < 4 fm), and the corresponding widths agree with experiment. The two states 
on the right are inconsistent with the 8i a = n/2 assumption since the two types of overlap calculations disagree at r < 4 fm. 
Computed widths of these two states also disagree with experiment. 



177 keV and 364 keV, respectively. 

There exist preliminary GFMC calculations of the mixing matrix element of the Hamiltonian for each of these 
doublets [591 . and this matrix element can be combined with the splitting of the doublet to determine its mixing 
angle (e.g. [82| )■ The splittings of the doublets are small enough that they are not resolved by the existing GFMC 
calculations, so the experimental splittings must be used to compute mixing angles. Since the recommended [68j | 165 
keV splitting of the 3 + doublet is less than twice the GFMC mixing energy, these two numbers cannot be combined 
to yield a real mixing angle. (Adopting a splitting of 310 keV [HJ produces good agreement with the experimental 
widths.) In the 1 + doublet, combining the GFMC mixing matrix element with the experimental splitting produces a 
poor match to the data (widths of 14.2 and 105 keV versus 10.7 ± 0.5 and 138 ± 6 keV from experiment). 

For each doublet, a mixing angle can be computed by % 2 minimization of the difference between theoretical and 
experimental widths. This produces good agreement with the measured widths, which was not guaranteed. The 
resulting mixing angle for the 3 + doublet has sin# m i x = 0.28, which is small relative to literature values 83]. The 
mixing angle from minimizing \ 2 for the 1 + doublet has sin(9 m ; x = 0.068, much smaller than both the value of 0.20 
from the GFMC mixing energy and the literature value of 0.21 [84|. 

Thus, attempts to compute separate widths for the upper and lower states in each doublet run aground on the 
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problem of finding mixing angles consistent with all available information. For each doublet, the quantity most 
independent of the mixing angle is the sum of the two widths. The sum is less sensitive to # m i x than are the 
underlying pole residues, but it is not quite independent of 6 m ix because the doublet states are not degenerate. It is 
the sum for each doublet that is shown in Table U and Fig. [31 with the T = 1 width computed in each case using the 
lower energy and the T = width the higher. These sums are in good agreement with experiment. 

C. Applications to recent measurements 

An important use of theoretical estimates of widths is in the identification of observed states. Most of the states 
considered here {A < 9, T < 1 MeV, dominated by nucleon decays) were found experimentally, and their spins and 
parities identified, long ago. Some exceptions are the first + and second 2 + states of 8 B and the entire low-lying 
spectra of 9 He and 9 Li. Here I attempt to shed light on these systems by calculating widths from VMC wave functions. 

1. 8 B 

Several theoretical models predict low- lying states of 8 B that have not been observed, as discussed in Ref. [7l|. It 
is possible that the states are simply too broad to be seen easily in experiments, and evidence was recently found for 
a + state at 1.9 MeV above the 8 B ground state [7l|. An i?-matrix fit to both elastic and inelastic scattering of 7 Be 
on protons indicated a + state with partial width 0.28 ± 0.14 MeV for decay to the 7 Be ground state and 0.33 ± 0.18 
MeV to the first excited state. Since these widths are within 2a of zero, I have not shown them in Table [U This state 
has some support from the calculation of Ref. [Tij , where it was found in computed phase shifts at the same energy 
using the merged no-core shell model and resonating group method. That calculation also indicates reason for caution 
in applying the integral relation: the phase shift does not approach ir/2, and indeed the VMC + state fails the ir/2 
consistency test discussed above. 

I show in Table U the predicted partial widths of the + state at its energies for the AV18+UIX and AV18+IL7 
Hamiltonians. Fig. [5] shows the dependence of this prediction on the assumed resonance energy. The claimed exper- 
imental widths are consistent with my results, though inconsistency of the overlaps with 7r/2 phase shift makes the 
significance of this consistency doubtful. I also include in Table Q] the AV18+IL7 energies and VMC-computed partial 
widths of the unobserved isobaric-analogue state in 8 Li. The 8 Li and 7 Li VMC wave functions are explicitly isospin 
rotations of the 8 B and 7 Be wave functions, and the 8 Li(0 + ) decays also fail the ir/2 condition on the overlaps. 

The authors of Ref. [7l| also find a 2 + state at 2.55 MeV excitation. Because this result has not been confirmed, 
I omit this state from goodness-of-fit statistics, but I have computed its partial widths to the 7 Be ground and first- 
excited states and included the results in Table fl] The overlap functions for these states are compatible with ir/2 
phase shift. However, the computed partial widths do not match those claimed in Ref. [71] : I find 0.51 ± 0.02 MeV 
to the 7 Be ground state versus 0.12 ± 0.04 MeV measured, and 0.039 ± 0.002 MeV to the 7 Be excited state versus 
0.24 ± 0.11 MeV measured. The origin of these differences is not clear. 

2. 9 He 

The spectroscopy of 9 He remains unclear despite several experimental studies |85T - [92l |. The ground state was 
originally thought to be a 1/2 - resonance state [85|. Subsequently, strong s-wave n- 8 He interaction was seen near 
threshold and argued to reflect a 1/2+ virtual stat e l8(j l87j . More recent experiments have found a smaller scattering 
length and thus less support for a virtual state [9~i1 . |92| . Because s-wave neutrons do not have true resonances passing 
through 7r/2 phase shift, I do not present a width for this state. 

Several other observations of resonances within a few MeV of the 8 He+n "threshold" have been claimed. (See 
Ref. [92( 1 for a summary.) Spin and parity assignments for all of these states are uncertain, and matching them to 
theoretical expectations has proven difficult [68j, |87|, 89] . Width estimates based on ab initio calculations could provide 
useful guidance, so I explore this possibility here. 

There have been four claims of a state around 1.2 MeV above the 8 He+n "threshold" [HI, |8(| [9^, [93| , and it has 
additional support from a study of possible analogue states [Hj]. This is generally assumed to be the lowest- lying 
p-shell state, with J T = l/2~, and there is conflicting information concerning its width. Ref. [85] found it to be 
narrower than their 0.42 MeV resolution; other experiments found 1 MeV 93], 0.10 ± 0.06 MeV [86j, and 2 MeV 
[9lj . As pointed out particularly by Barker [8t|, it is difficult to reconcile widths considerably narrower than 1 MeV 
with the expected strong single-particle character of the l/2~ resonance. Theoretical calculations also place the 
l/2~ state a few MeV higher (e.g. at 3 and 4 MeV in GFMC calculations with Illinois-6 and Illinois-2 three-body 
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FIG. 6: (Color online) Predicted partial widths as functions of assumed resonance energy (relative to the 7 Be ground state 
threshold) for the first + state in 8 B. The thicknesses of the bands indicate the errors from Monte Carlo sampling. The band 
bounded by solid curves indicates the partial width to the ground state of 7 Be, and that bounded by dashed curves shows the 
partial width to the first excited state. The boxes (offset from each other slightly in the horizontal direction for visibility) show 
the experimental results of Ref. [7lJ without a correction from J?-matrix formal width to observed width. (This can amount 
to a ~ 30% reduction at the high end of the allowed width.) Also indicated is the best estimate of the resonance energy from 
GFMC calculations using the AV18+IL7 Hamiltonian. 

forces, respectively). One possibility is reduction of the l/2~ energy by an sd-shell component that is missing in the 
calculations pj|. 

Additional resonances have been found at higher energies, one around 2.3 MeV [85|, HH, [92| with claimed width 
0.7 ± 0.2 MeV, one around 4 MeV [8rj[fl|. and another around 5 MeV [H HI]. In addition to the 1/2" state, a 
3/2~ state is expected theoretically, though also at higher energy. It is also likely that sd-shell intruder states with 
J v = 5/2+ and 3/2+ are present in the low- lying spectrum. 

I computed widths of p-shell states with = l/2~ and 3/2~, but the results proved difficult to match with 
experimental data. The short-range overlaps of the 1/2 - state computed directly and by the integral method are 
in nice agreement. Since it is unclear what resonance energy should be used in the calculation, I have computed 
widths from the VMC wave function using a range of channel energies in the integral relation. Fig. [7] shows these 
results. They mainly demonstrate the argument of Barker [89| that there is a mismatch between the narrow width 
of the l/2~ state claimed by Bohlen et al. [86| and theoretical expectations of a strongly single-particle state. The 
resonances claimed by Belozerov et al. [Hf and Golovkov et al. [9l| are consistent with the computed l/2~ width, 
but that assignment makes interpretation of the 1.3 and 2.4 MeV states of Bohlen et al. [86| difficult: at least one 
would have to be an sd-shell state. 

All channel energies below 6 MeV for the 3/2~ state yield computed widths of less than 5 keV. This could match 
either of the states at 4.3 and 5.3 MeV in Ref. (86[, which were found to be narrower than the 100-keV experimental 
resolution, but the VMC overlaps of this state are inconsistent with 7r/2 phase shift. For this reason, the integral- 
relation results are probably not reliable predictions of the width. 

VMC wave functions also exist for states with J v = l/2 + ,3/2 + , and 5/2 + , but these wave functions with sd-shell 
components have not reached the same level of development as the VMC p-shell states and have not been published. 
As mentioned above, the 1/2+ state should not be observable as a resonance. The VMC wave functions for the 
3/2+ and 5/2 + states indicate that they are mainly made by coupling the 2 + state of 8 He to an s-wave neutron and 
therefore also not true 6 a — tt/2 resonances. I did not attempt calculations of partial widths for decay to the 8 He(2 + ) 
state. Calculations of the partial widths for decay from positive-parity states to the 8 He ground state by emission of 
a d-wave neutron yielded partial widths of less than 60 keV and overlaps inconsistent with 7r/2 phase shift. 



21 




1 2 3 4 5 6 

E R (MeV) 



FIG. 7: (Color online) Boxes (with sizes indicating reported errors) show claimed He resonances that have reported widths. 
The band labeled "VMC" shows the predicted width of the l/2~ state as a function of the resonance energy assumed in the 
integral relation, and its width reflects the statistical error of the Monte Carlo integration. Where no error was reported for 
a width (for example, Ref. [9ll ] has "~ 2 MeV" for the width of the state at 2.0 MeV), I have assigned an error of 0.5 MeV 
and indicated this with broken lines at the upper and lower limits. I have omitted some very broad states with missing or 
lower-limit errors: a state with E R ~ 3 MeV and V ~ 3 MeV from Ref. [93| and one with E R > 4.2 MeV and F > 0.5 MeV 
from Ref. [9lJ. See Ref. [92| for additional reports of states without measured widths. The displayed widths are from Seth 
et al. [85|, Belozerov et al. |93j|. Bohlen et al. [8(|, and Golovkov et al. 91]. Downward arrows indicate widths that include 
~ 400 keV instrumental resolution. 

3. 9 Li 

The most recent data compilation for A = 9 lists five low-lying states of 9 Li and only assigns firm spins and parities 
to the lowest two [68J. A more recent paper [63[ identifies the third state as J w = 5/2 ~ by comparing spectroscopic 
factors of VMC wave functions (older versions of those used here) with measured (d, p) cross sections. Those authors 
also assigned J 71 = 3/2~ and 7/2~, respectively, to the next two states. This assignment was based partly on the 
ordering of states in theoretical calculations and partly on the assumption that widths should correlate with computed 
spectroscopic factors. The results presented in Tableland Fig. [3] support these assignments by reproducing the widths 
of all three unbound states. 

D. Comparison with other width estimates 

I conclude by considering other ways to estimate widths from VMC wave functions and comparing them with the 
integral method. In the absence of integral relations or explicit scattering calculations, widths must be estimated 
from spectroscopic factors. In applications of the shell model, one often assumes that the width is the product of the 
spectroscopic factor, 

/>oo 

S a = / [R a (r)f r 2 dr , (53) 
Jo 

and the single-particle width. VMC spectroscopic factors might in fact be more suited to this procedure than those 
from a shell model, because shell models are typically confined to a single value of the principal quantum number, 
while the large amount of correlation in VMC wave functions guarantees that Eq. (|18p picks up contributions from 
all major shells. 

The crudest estimate of the single-particle width is the "Wigner limit" [§4j]. On the basis of a causality argument, 
the width of a resonance can be shown to have an approximate upper limit of 

T w =2-^P z (fc a fc), (54) 
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FIG. 8: (Color online) Comparison between experimental widths and widths estimated as the VMC spectroscopic factor times 
the Wigner limit. The vertical axis shows the ratio of these numbers, using a radius of 4 fm to compute the Wigner limit. 
Symbols are as in Fig. O and the labels from that figure may be used to identify states here. Comparison with Fig. [3] indicates 
that the integral relation is a significantly better predictor of widths than S a F\y- Note the different vertical scale from Fig. 



where b is a "matching radius" (typically ~ 4 fm) defining the edge of the interaction region, and Pi is the penetration 
factor of Ref. |40| . Since Ty/ is (approximately) an upper limit on the width that a resonance can have, it might 
approximate a single-particle width. Some authors define T\y to include an additional numerical factor multipl ying 



Eq. (|54|) , reflecting assumptions about the wave function inside the interaction region. Teichmann and Wigner '94 1 
assumed a constant wave function to arrive at a factor of 3/2. Other authors make more elaborate assumptions and 
arrive at a factor of (21 - 1)/(2Z + 1) for I ^ 95]. I take Eq. flM]) to define Tw- 
in Fig. [8] I use the Wigner limit to estimate widths of the states under consideration. In each case, I multiply 
for a 4 f m radius by the VMC spectroscopic factor from Eqs. (fP5)) and ((551) . These estimates plainly do not reproduce 
measured widths as well as the integral relation. The weighted mean ratio of "theoretical" to experimental width for 
this method, restricted to states consistent with 7r/2 phase shift and purely nucleon-emission decay, is 2.49 ± 0.52; 
the reduced x 2 is 1845 for eight degrees of freedom. (Recall that (ri nt egr a i/r e xpt) = 1-08 ± 0.04 with xl = 1-6 for 
the same set of states.) The mismatch between S a Tw and r expt can be reduced by choosing a smaller numerical 
factor to define the Wigner limit, but that does not remove the large scatter in S a Tw /r cX pt- Better agreement with 
experiment is likewise achieved with a smaller radius, but good agreement requires an unphysically small radius in 
the neighborhood of 2 fm, which again does not remove the large scatter in the ratio S a Tw /r eX p- Estimates from 
Tw are typically ~ 20% smaller if they are estimated as "observed widths" from Eq. (|5ip , using 

h 2 

la = Sa-j- • (55) 

This helps significantly with neither the overall scale nor the large scatter of the predicted widths. 

A better estimate of the single-particle width, and one perhaps more widespread in shell-model studies, is based 
on Woods-Saxon potentials. One assumes a potential well of "standard" radius and diffuseness and adjusts its depth 
to produce a resonance at the correct energy. The width Tws of this resonance is then taken as an estimate of 
the single-particle width. A range of geometric parameters for the potential is usually considered, because the most 
appropriate values are not known a priori. 

I computed estimates of this kind, using diffuseness 0.65 fm and Woods-Saxon radius 1.1 fm x (A — l) 1 / 3 , with A 
the mass number of the resonance state (the defaults in a code provided by B. A. Brown [96]). I neglected variation 
of these parameters and estimated the single-particle width T-ws from the FWHM of the peak in the computed cross 
section for the given Woods-Saxon well. The comparison of SaFws with r cxpt is shown in Fig. [S] The Woods- 
Saxon estimates are systematically low, with a weighted mean (S a T-^s/^expt) = 0.74 ± 0.15 for the same eight cases 
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FIG. 9: (Color online) Widths estimated as the VMC spectroscopic factor times the Woods-Saxon single-particle width Tws, 
divided by experimental widths. Symbols are as in Fig. [3] and the labels from that figure may be used to identify states here. 
Comparison with Figs.[3]and[8]indicates that this approach predicts experimental widths more accurately than the Wigner-limit 
approach but still not as well as the integral relation. 



considered above and xt = 34 for eight degrees of freedom. I conclude from this exercise and the similar exercise 
using Wigner limits that VMC widths computed by the integral method are more useful predictors of experimental 
widths than are the VMC spectroscopic factors. 



V. SUMMARY 



I have presented plausibility arguments, supported by detailed derivations in the literature, that widths of resonant 
states can be estimated by evaluating an integral over pseudobound ab initio wave functions. This approach is 
approximate, but it avoids a great deal of computation and human labor that would be needed in full-on scattering 
calculations and would often be complicated by coupled channels. It is nicely suited to quantum Monte Carlo 
calculations in that it is insensitive to the difficult-to-compute tails of the many-body wave functions, it involves a 
short-range integral amenable to Monte Carlo integration, it uses more information about the Hamiltonian than is 
encoded in the variational wave function, and it can be applied to resonances narrower than the practical energy 
resolution of the GFMC technique. Related integrals yield overlap functions for bound states, and these overlaps are 
guaranteed to have the correct shapes in their long-range asymptotics even when the variational wave function does 
not. These may be useful for calculations of spectroscopic factors and of transfer and knockout cross sections. 

I have implemented integral-method width calculations for one-nucleon emission from wave functions computed by 
the variational Monte Carlo method. It yields widths in good agreement with experiment for several states in the 
7 < A < 9 mass range. Cases of disagreement always involve either open channels for which I have not accounted or a 
resonant wave function that is not strongly peaked in the interaction region. I have shown that widths predicted in this 
way are closer matches to experiment than are naive combinations of ab initio spectroscopic factors with Wigner-limit 
or Woods-Saxon estimates of single-particle widths. The integral method is thus a useful tool for estimating widths 
from ab initio methods that produce pseudobound wave functions. 

For the longer term, the calculations presented here represent a learning problem for application of integrals of 
the type in Eqs. (|3"2")) - (j3"3"]) to QMC wave functions. Application to GFMC wave functions of the methods used here 
will be straightforward and mainly involve additional bookkeeping similar to that used in Ref. [6fj| for direct overlap 
calculations. Integrals of the kind considered here are likely to find their most extensive use in ab initio calculations 
of coupled-channel scattering and reactions, and a major goal of the work presented here is to prepare the way for 
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such calculations. 
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